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Abstract. A sequence of completely positive maps can be decomposed into 
quantum trajectories. The geometric phase or holonomy of such a trajectory is 
delineated. For nonpure initial states, it is shown that well-defined holonomics 
can be assigned by using Uhlmann's concept of parallel transport along the 
individual trajectories. We put forward an experimental realization of the 
geometric phase for a quantum trajectory in interferometry. We argue that 
the average over the phase factors for all quantum trajectories that build up a 
given open system evolution, fails to reflect the geometry of the open system 
evolution itself. 
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1. Introduction 

Imagine a quantum particle prepared in a pure state. If the state evolves in a 
unitary fashion so that it remains pure, this particle is said to constitute a closed 
system. Berry [^| and others [El El El have shown that a closed system picks up a 
phase factor of geometric origin; a geometric phase factor, or quantum holonomy. 

But what happens if the system fails to be closed? Is there a natural gen- 
eralization of the geometric phase to open systems? Considerations of phases of 
geometric origin for open quantum systems go back to work by Uhlmann [ |S] on 
holonomy accompanying density operators. Garrison and Wright have assigned 
complex- valued geometric phases to open systems on a phenomenological basis by 
the use of non-Hermitian Hamiltonians. Ellinas et al. [[7] have introduced geometric 
phases associated with adiabatic evolution of eigenmatrices of the Liouvillian su- 
peroperator describing an open system. Gamliel and Freed demonstrated that 
the adiabatic Berry phases may appear also in the presence of open system effects, 
as phase factors in the off-diagonal elements of the density operator expressed in 
the instantaneous energy eigenbasis (see also [|H1 llUp. 
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More recently, open system geometric phases have been revisited from various 
perspectives, such as interferometry 1121 II quantum trajectories H5j . 
phase distributions [E|> the adiabatic theorem and stochastic Schrodinger 

evolutions [I18|. A reason for this renewed interest is the need to better understand 
the robustness properties of geometric phases in the presence of open system ef- 
fects, triggered by the conjectured importance of quantum holonomy for fault 
tolerant quantum computation. To address the robustness issue, it is commonly 
assumed that the system evolves under a Markovian master equation. The quan- 
tum trajectory approach is specifically designed for the analysis of the robustness 
of geometric phases to Markovian open system effects (see, e.g., [ 1201 121|). 

In this paper, we focus on the quantum trajectory approach to the open sys- 
tem geometric phase. The purpose is fourfold. First, we wish to reformulate the 
quantum trajectory approach to the geometric phase on the basis of discrete sets of 
completely positive maps (CPMs), without invoking continuous master equations. 
This reformulation is a generalization of previous treatments [^J El ( see also [I22p 
in the sense that it yields the master equation in the continuous limit. Secondly, 
we wish to extend [ 114] to quantum trajectories that start in a nonpure state. To 
this end, it seems most natural to use Uhlmann's concept of parallel transport [E] 
along the trajectories. Thirdly, we put forward an experimental procedure to imple- 
ment the geometric phase for an individual trajectory. Finally, we briefly address 
whether the geometric phases for individual quantum trajectories can be used to 
define a meaningful geometric phase for an open system. The paper ends with the 
conclusions. 

2. Open quantum systems and quantum trajectories 

Consider a quantum system s in contact with some environment e, also of quantum 
nature. Suppose s + e is in the state g. Any observable property pertaining only 
to s is given by the reduced density operator p obtained by tracing g over e. Let 
U(t, to) be the time evolution operator for s + e. Further, let |eo), ■ ■ ■ , le^) (/x could 
be finite or infinite) be a complete basis for e and suppose that g = |eo) (eo| <8> p is 
prepared at to- Then, for an arbitrary t, the state of s reads 



where E p (t,to) = (e p \U(t, to)|eo) are the Kraus operators [ [231 that constitute a 
Kraus representation of £t,t - The map £tt 1S trace preserving and completely 
positive, i.e., it takes normalized density operators into normalized density operators 
and all trivial extensions 1 ® £t,t a likewise. 

To introduce the concept of quantum trajectories it is convenient to consider 
N identical copies of the environment. Prepare them in the product state |eo) = 
|eg) <8> . . . <S> le^" 1 ), with the superscript labeling the copies. Let s interact only 
with the first copy between £q and t\, only with the second copy between t\ and t2, 




(1) 
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and so on. These interactions are described by the unitary operators U°(t, to), . . . , 
C/ Ar_1 (t, fiv-i), respectively. It results in the composite map 

p->£(p)= £ tN ,t N -x o . . . o S tl , to (p). (2) 

Here, E p (tk+i,tk) = (e£|[7*(£fc+i, £fc)l e o)> P = 0,...,p, can be taken as Kraus 
operators for £t k+1 ,t k , k = 0, . . . , N — 1. In terms of the above s + e picture, the 
relevant part of the input state at tk is |e§) (e§ | <g> pk, where po = p and Pke[l,N-l] = 
£t k ,t k _ 1 ° • ■ ■ £ti,t (p)- ^ quantum trajectory a is a sequence of (unnormalized) 
states p — ► pf — > . . . — > p^v, where 

Pfe - ^(fe) (*fc, *fc-i) ■ • • E a{1) {tut^pE^ (t!, t ) . . . £7+ (fe) (t fc , tft-i), (3) 

for fc = 1, . . . , N . Here, a (I) € [0, ^t] is the li/i element of a sequence of indexes. The 
map £ is recovered by summing over all trajectories. 

3. Geometric phase of quantum trajectories 

There is in general an infinite number of possible Kraus representations for a given 
CPM. For the composite map £ in Eq. J2J, different choices of Kraus representations 
in each step lead to different sets of trajectories. It follows that no physical meaning 
can be associated to an individual trajectory without imposing some additional 
physical constraint. 

Nevertheless, it is still possible to formally associate a geometric phase or holon- 
omy to an individual quantum trajectory. Here, we wish to do so. In particular, we 
wish to stress the importance of whether the initial state of the trajectory is pure 
or nonpure. In the pure case, which has been treated in [ 1141 122) . the standard Pan- 
charatnam connection can be used, while in the nonpure case we need to address 
the issue of parallel transport of density operators. 

Let us first consider the case where the initial state is pure, i.e., p = |V')(V , I- 
Here, the trajectories may be lifted to sequences of (unnormalized) vectors: \ip) — > 
\ip%), where |^) = E a{k) (t k , t fc _ x ) . . . E a{1) (h, t Q )\^). The Pan- 
charatnam connection [ |2] yields the geometric phase factor of the trajectory a 
as 

7° = $[{^>{^-i>."WW, (4) 

where &[z\ = z/\z\ for any nonzero complex number z. 

When the initial state is nonpure, though, one needs to consider geomet- 
ric phases for mixed states, since the trajectories then contain nonpure states 
p,pf,... , p% (see Eq. ©)■ To deal with this, it seems most natural to use the 
Uhlmann approach [5, since this approach is particularly well adopted to dis- 
crete sequences of density operators. Assume that p,pf,... , p^ are full rank and 
introduce the amplitudes W — \fp~V and W£ — \fitftV k , k = 1,...,^, where 
V, V", ■ ■ ■ , V§ are the unitary operators or 'phases' of the amplitudes. For a given 
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but arbitrary V, V" is fixed by the parallelity condition W^W > and similarly 
V^ +1 is fixed iteratively by W^^W^ > 0, for k = 1, . . . , N— 1. In this way, the final 
unitary operator V$ is uniquely determined by the trajectory a up to the arbitrary 
phase V of the initial amplitude W. Explicitly, 

x( v / pfp?v / pff ^ ''^^(Vpfp^pf) 11 ' 'Vfiy/pV. (5) 

To remove the arbitrary phase V, define 

U a =V£V\ (6) 

which is the holonomy of the trajectory a. Note that U a reduces to 7" in the limit 
of pure initial states. 

4. Accessing the geometric phase of a quantum trajectory 

As already mentioned, no physical meaning can be associated with a single tra- 
jectory and its concomitant geometric phase or holonomy if no further physical 
constraint is imposed. This is so because there are infinitely many Kraus represen- 
tations of a given CPM, and therefore infinitely many equivalent ways to decompose 
the open system evolution into quantum trajectories. Nevertheless, the geometric 
phase for a single quantum trajectory can be accessed in principle by performing 
measurements on the environment. 

To sec this, let us first make some simplifying assumptions. Suppose that 
p = IV'XV'I (pure initial state) and that the observables O^, k = 0, . . . ,N — 1, with 
eigenstates |e§), . . . , \ez), are measured projectively. We further assume that these 
measurements are performed precisely at t\, . . . , t/v, i.e., that the first copy of the 
environment is measured at ti, the second at t 2} and so on. In the first step, s + e 
evolves as (ignoring all copies except the first one) 

leg) (e° 1 ®|V> M ^ [/°(ti,to)|eg)(eg|®|V)^|C/ 0t (ti,io) 

- |e°(i)}(e° (1) |®^ Q (i)(ii,io)|V)^l^ (1 )(*i^o) (7) 

for the registered outcome a(l) G [0,/i]. In the second step, the state of s + e 
similarly evolves as (now ignoring all copies except the second one) 

(el I ® E a{1) (h ,toM) (VK (1) (*i , to) 
-> V\t 2l t x )\e\) (ej| ® E a{1) (h, i )|V) (^1^(1) *o)tf lf (*2, *i) 

- |ea(2))(e«(2)l®^a(2)(i2,tl)^a(l)(il,io)|V')^l^ (1 )(tl,to)4(2)(*2,il) (8) 

for the registered outcome a(2) S [0, yu]. Continuing in this way provides a physical 
realization of the trajectory a. 
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We now show that this measurement procedure can be used to experimentally 
implement 7°. Suppose is the input state of the internal degree of free- 

dom (e.g. spin) of an ensemble of particles in an ordinary two-beam Mach-Zchndcr 
interferometer. The evolution in Eq. Q is implemented in one of the interferom- 
eter arms by post-selecting those particles for which a(l) is registered. A variable 
U(l) phase shift is applied to the other arm. The maximum of the resulting in- 
terference fringes is obtained for $ [(V'l-E^i) (ti, ^o)|^)] — [(V'i" l^)] ■ In the next 
step, the internal state E a ^ (t\, to)\ifi) (V'l-^am (ti, to) is used as input to a second 
interferometer in which Eq. © is implemented in one arm by post-selecting the 
outcome a(2). The other arm is again exposed to a variable U(l) shift, yield- 
ing an interference maximum at &[(ip\E^ a ^ (t±, to)E^ 2 ^ (ta, ti)E a ^ (ti, to)|^)] = 
$[(^2 l^i)] • Continuing in this way up to ip^ and back to ip yields the phase 
factors $[(V#^)],$[(V$|V$_i)]> ■ ■ - ,$[(^1 1^)]- By taking the product of these 
phase factors, we obtain 

*[<^l^>]^[<^l^-i>] ---^[WIV')] = ^[<Vi^X^I^-i>-.-<^?l^>], (9) 

which is precisely the geometric phase factor j a in Eq. J3J. It should be noted, 
though, that in order for this procedure to work it is in each step necessary to erase 
any path information related to the state change of the environment caused by the 
measurement. 

5. Open system geometric phase 

So far, we have discussed how a geometric phase factor 7" or holonomy U a can be 
assigned to a single quantum trajectory a and how 7" can be verified in interferom- 
etry. However, since an individual quantum trajectory can differ very much from 
the original evolution £, it is not obvious in what way 7° gives information about 
the geometry of the path that is generated by £ itself. Such a connection can be 
established if it is possible to use the geometric phases for the trajectories to assign 
a well-defined geometric phase for £. Here, we briefly address this issue. 

The map £ is recovered by summing over all trajectories a. In the case of pure 
initial state \tp), it is therefore natural to guess that there is a geometric phase 
factor of £ that is the average over the phase factors of the individual quantum 
trajectories. Explicitly, 

T £ =Y J P a l a , (10) 

a 

where p a — (iPnI^n) ^ s the probability of a. Unfortunately, though, this quantity 
depends on the specific decomposition of £ into trajectories. Thus, unless a specific 
set of trajectories is singled out (e.g., by repeatedly measuring the environment, 
as described in the preceding section), no physical meaning can be assigned to Tg. 
An analogous conclusion in the context of stochastic Schodinger equations has been 
drawn in [ I18j . 
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6. Conclusions 

The quantum trajectory approach has recently been proposed to tackle the issue 
of geometric phase [ 1141 I22| and quantum holonomy [ I15| for an open quantum 
system. The main merit of the approach appears to be its usefulness in the analysis 
of the robustness of geometric phase based quantum gates to open system effects 
[ 1201 121) . Here, we have reformulated the quantum trajectory approach to the 
geometric phase, on the basis of sequences of completely positive maps. 

It has been stressed [ 1141 [22] that the quantum trajectory approach is useful in 
that it avoids the apparent need to introduce parallel transport of nonpure states. 
This holds only under the assumption that the trajectory starts in a pure state. 
We have pointed out that even in the case where the initial state is nonpure it is 
possible to assign well-defined holonomies by using Uhlmann's concept of parallel 
transport [ 5 along the quantum trajectories. 

The geometric phase of an individual trajectory can be made physical by re- 
peatedly measuring the environment of the system. This opens up the possibility to 
experimentally implement such phases by measuring post-selectively the environ- 
ment in one of the arms of a standard Mach-Zehnder interferometry. In this way, 
we have put forward an iterative procedure to measure the geometric phase of a 
single quantum trajectory. 

Finally, we have addressed the question whether the geometric phases for the 
trajectories comprising an open system evolution can be used to define the geometric 
phase of the open system evolution itself. A natural approach is to sum over the 
individual phase factors weighted by the corresponding probabilities. However, this 
works only if a specific set of trajectories is singled out by some physical means. 
It is interesting to note the analogous situation for the geometric phase proposed 
in [E]- This phase depends on the additional detailed knowledge of the system- 
environment interaction. Thus, neither the approach in [ lllj nor the averaging over 
quantum trajectories put forward here, result in a genuine concept of geometric 
phase for an open system. It seems more appropriate to use the path of density 
operator itself as the basis for such a concept [03 I13| . 
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